Role of Hidden Slow Degrees of Freedom in the Fluctuation Theorem 
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The validity of the fluctuation theorem for entropy production as deduced from the observation 
of trajectories implicitly requires that all slow degrees of freedom are accessible. We experimentally 
investigate the role of hidden slow degrees of freedom in a system of two magnetically coupled 
driven colloidal particles. The apparent entropy production based on the observation of just one 
particle obeys a fluctuation theorem-like symmetry with a slope of 1 in the short time limit. For 
longer times, we find a constant slope, but different from 1. We present theoretical arguments for 
a generic linear behavior both for small and large apparent entropy production but not necessarily 
throughout. By fine-tuning experimental parameters, such an intermediate nonlinear behavior can 
indeed be recovered in our system as well. 
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Introduction. — Basic concepts of thermodynamics and 
statistical physics implicitly rest on a separation of all 
degrees of freedom into observable and non-observable 
ones. Heat exchange, for example, is associated with the 
myriads of fast degrees of freedom which are not resolved 
dynamically, whereas work typically involves a few con- 
trolled, slow degrees of freedom. In systems without a 
clear-cut time-scale separation, ambiguities and inconsis- 
tencies may arise if such concepts are still applied naively. 
Here we explore this issue for one of the arguably most 
relevant concepts, entropy production Astot, for non- 
equilibrium steady states (NESS). For such states, the 
fluctuation theorem (FT) refers to a remarkable symme- 
try that quantifies the probability p of observing trajec- 
tories with negative total entropy production as 



ln[p(Astot)/p(-Astot)] = aAsto 



with a~l and Boltzmann's constant set to unity jlH5|- 
The FT has been proven for two types of dynam- 
ics. First, for deterministic dynamics the proof rests on 
the chaotic hypothesis, including time reversibility and a 
phase-space contraction associated with dissipation 
Second, in stochastic dynamics the FT requires the con- 
cept of entropy production along trajectories and can be 
proven for Markovian systems [3|-|5| . The latter dynamics 
applies to experiments on driven colloidal particles [SUtI 
and a harmonic oscillator coupled to a thermal bath [81 . 
Experimental tests of FT-like symmetries have also been 
reported for Rayleigh-Benard convection turbulent 
flow [lol , granular matter [ll| , and self-propelled parti- 
cles [1^. For these systems the appropriate class of dy- 
namics is less obvious and hence the status regarding the 
assumptions of the FT is unclear a priori. One should 
also appreciate that the measured observable for some of 
these systems is typically not Astot directly, but rather 
some dimensionful quantity, like, e.g., the injected or dis- 
sipated work [iil, m, which requires a temperature for 
a unique conversion to entropy. Strictly speaking, the 



FT is thus verified only if this temperature can be deter- 
mined independently and if it leads to a = 1. Comple- 
mentary, in more recent reports, the validity of the FT 
(with a = 1) is assumed and used to gain information 
about such a dimensionful factor connecting the actual 



observable with entropy production 

For stochastic dynamics, the proof of the FT with 
a = 1 rests on a time-scale separation. Fast degrees of 
freedom contribute to an effectively white noise leading 
then to a Markovian dynamics of the slow degrees of free- 
dom. Entropy production can be deduced from observing 
the dynamics of all slow degrees of freedom. If some of 
these degrees of freedom are not, or cannot, be observed 
the inferred entropy production is only an apparent one 
for which the status of an FT-like symmetry is unclear a 
priori. Theoretical efforts to describe coarse-graining in 
(1) general have been restricted so far to the case of well sep- 



arated time scales 



16 - 1^ , and how such coarse-graining 
affects bounds on dissipated work [l^ . In the framework 
of electronic devices FT-like symmetries for currents have 
been discussed in Refs. 

In the present paper, we investigate the role of hidden 
slow degrees of freedom on apparent entropy production 
for a paradigmatic system with two magnetically coupled 
driven colloidal particles. 

Apparent Entropy Production. — The total entropy pro- 
duction Astot is given by the sum of the entropy changes 
of the heat bath and system 



Astot = + hi [p^(xO)/p^(x*)] 



(2) 



where Q = /g cJt ^"^j^ a;i(r)Fi(x(r)) is the heat transfer 
of all n degrees of freedom x = (x'l . . . x„) to the solvent 
at temperature T . Here, Xi is the actual velocity and 
Fi is the total force acting on the z*^ degree of freedom. 
The change of the system's entropy includes the station- 
ary probability distribution p'^(x°) (p'^(x*)) of finding the 
initial (final) state of the system along the trajectory of 
length t. In a NESS, the system satisfies the station- 
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ary Smoluchowski equation = —J^i^^xiJii^), where 
the probabiHty current j^x) = p^{x)L'f{x.) isgiven as 
the product of p'^ and the mean local velocity Q 



z.f(x) = i^o m^)/T~d^^ Inp^(x) 



(3) 



with Do the bare difFusivity. Multiplying Eq. ([3]) with 
Xi and integrating over time yields the total entropy pro- 
duction as given in Eq. ([2]), 



Astot = J dTY,it{'rXi^ir))/Do, 



(4) 



where the sum involves all n degrees of freedom x. If 
only the first h of these, x = {xi . . .Xn), are accessible, 
an observer is forced to deduce all information from these 
trajectories. The actual velocities x can still be measured 
correctly, whereas the mean local velocity obtained from 
the accessible trajectories is 



i',^(x)p^(x|x)(ix 



(5) 



the conditional probability p''(x|x) for x = 
Xn) at fixed X. Hence, the apparent entropy 



with 

(Xfi+l, 

production becomes 



dT^i,(r)i>ni(r))/i?o, (6) 
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FIG. 1: (color online) (a) Schematic representation of the 
system and measured tilted potentials (Vi" = V2 = 181 ksT 
and /i = -/2 = 57 ksT/^m calculated via (b)-(d) Sta- 

tionary probability distribution p^{xi,X2) for different plasma 
parameters F. The white circles indicate typical trajectories. 



where the sum runs over the h accessible degrees of free- 
dom only. In this Letter, we investigate the conditions 
under which this quantity obeys a FT-like symmetry. 

Experiment. — We have created two non-overlapping 
toroidal traps with radius R = 3.5 iJ.m and a center- 
center distance of 17 |J.m by a single laser beam (A = 
1070 nm) which was deflected on a galvanometric mirror 
unit (for details refer to 0, [ISl)- Each trap contained 
a single paramagnetic colloidal particle with a 2.6 |i.m 
radius (Microparticles, Berlin) and labeled by an index 
J = 1,2 (see Fig. 1(a)). The traps are approximately 
50 |i.m away from the lower surface; therefore hydrody- 
namic interactions with the walls are negligible 25| . The 
scanning frequency was adjusted to 41 Hz, which leads 
to quasistatic tangential forces fi acting on the i*^ par- 
ticle along the toroidal traps whose amplitude depends 
on the laser intensity. In our experiments, the time for a 
full revolution of each particle was adjusted to 10 s. Syn- 
chronized with the scanning motion of the laser beam, 
its intensity was sinusoidally modulated with an acousto- 
optic modulator which finally leads to an effective optical 
potential Ui{xi) = Vi{xi) — fiXj with Vi{xi) = sin(xi), 
where Xi is the particle position along the trap circum- 
ference in units of ttR (see Fig. 1(a)). Accordingly, both 
particles reach NESS, where fi and can be controlled 
independently. 



A coupling between the two NESS is obtained by a 
static homogeneous magnetic field B applied perpendic- 
ular to the sample plane. This field induces magnetic 
moments m to the particles leading to a repulsive dipolar 
particle interaction W{xi,X2) = (/io/47r) TO^/r"^(xi, 2:2). 
Here, ^0 is the magnetic constant and r the particle 
distance. For small magnetic fields {B ^ 40 mT), 
as in our experiments, the magnetic moment scales as 
TO « (7TOo/3)-B with 7 = 30 T~^ and toq = 5.9 x 
10"^"^ Am^ [28|. The strength of the dipolar coupling can 
be conveniently characterized by a dimensionless plasma 
parameter F = AW/(A:bT), where AW corresponds to 
the difference of the coupling at the smallest and largest 
particle distance r. Additional particle interactions, e.g., 
hydrodynamic coupling or optical binding, are negligi- 
ble at the chosen trap separations as confirmed by the 
independent motion of the particles in the absence of a 
magnetic field. 

Results. — We investigate the effect of coupling by 
preparing two identical NESS, where the potential min- 
ima face each other. Figures[ljb)-(d) show the stationary 
probability distribution p^{xi,X2) as color coded back- 
ground and as white circles one exemplary trajectory. 
In the absence of coupling, the peak in p^ corresponds 
to the flattest part in the potentials, where the parti- 
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FIG. 2: (color online) (a) Distribution of the apparent en- 
tropy production p(Astot) for different trajectory lengths t 
and plasma parameters F. (b) Section of previous histograms 
around Astot = 0. (c) Corresponding In [p(Astot)/p(-Astot)] 
as a function of Astot- The dashed black line has the theo- 
retically predicted slope 1, whereas the red line is a linear fit 
with slope Q = 0.65. 



clcs slow down and therefore are most likely to be found 
(see Fig. [Ijb)). Under strong coupling conditions (see 
Fig. [Ijd)), the repulsion hinders the particles from com- 
ing close to each other. The peak of vanishes since the 
approach of one particle kicks the other one away, lead- 
ing to motion like a Newton's cradle. In the intermediate 
regime (see Fig. [ijc)), the coupling interaction is com- 
parable with the energy loss acquired while the particle 
moves along half a circle. Here the full interplay between 
drift, diffusion, and interaction has to be taken into ac- 
count and we expect the influence of hidden degrees of 
freedom on the FT to be most prominent. Therefore, we 
concentrate on Astot associated with the motion of only 
the first particle, which represents the observed degree 
of freedom, whereas the coupling allows us to uniquely 
control the infiuence of the second particle, which acts as 
a hidden degree of freedom. 

The black histograms (closed bars and line) in Fig.[2l^a) 
show the distribution of the apparent entropy production 
p(Astot) in the absence of coupling obtained for trajec- 
tories of length t = 1.75 s and 10 s, respectively. The 
peaked distribution shifts with elapsing time to the right 
with peak height maxima occurring at positions which 
correspond to the energy loss associated with full revolu- 
tions of the particle, 27ri?/ = 1250 keT. To investigate 
the FT, rare events with negative entropy production 
have to be sampled with high accuracy. This constrains 




FIG. 3; (color online) (a) Slope a vs plasma parameter F 
for t — 1.75 s. (b) Slope a for different trajectory lengths 
t. The black squares correspond to F = and the red dots 
to F = 300. The deviation of the black squares from a = 1 
(black dashed line) determines the statistical errors to be less 
than 5%. 



the maximal trajectory length t to approximately 2 s and 
the range within the FT can be tested to ±3. Figure [IJb) 
shows this section of the black (closed bars) histogram 
around Astot = 0. The excellent agreement between the 
logarithm of the probability ratio p(Astot)/p(— Astot), 
black squares in Fig.^^c), and the black dashed line with 
a slope of 1 confirms the validity of Eq. ([1} for uncou- 
pled states. The red histogram (open bars) in Figs. HJa) 
and HJb) demonstrates the situation for coupled states. 
Most prominent is the enhanced probability at Astot = 0. 
Since the red dots in Fig. ^c) do not agree with the 
dashed line of slope 1 this apparent entropy production 
does not obey the FT. Rather a linear relation as given 
by Eq.([T|) with a ~ 0.65 is found. 

In additional experiments, we observe linear relations 
according to Eq. ([T|) with different slopes a, which de- 
pend on two parameters: (i) the plasma parameter F, 
and (ii) the trajectory length t, as shown in Figs, ^a.) 
and ^h). Clearly, the FT is confirmed for arbitrary 
trajectory lengths in uncoupled states (black squares in 
Fig. E^b)). The obvious dependence of a on F resembles 
the transition from an uncoupled to a coupled state. For 
F = 300, a decays with increasing length t, from 1 to 
0.65. A similar time dependence has been also observed 
in |14| : however, there it was not identified as an inherent 
feature of hidden degrees of freedom. In additional mea- 



surements performed for coupling two different NESS, we 
also found such a linear relation. Therefore, we exclude 
symmetry as the sole origin of this behavior. 

Discussion. — First, we explain why for t — >■ the slope 
a approaches 1. In general, deviations from the FT must 
be caused by the interaction with the hidden particle. In 
a short time expansion to lowest order in t, we can ne- 
glect changes in the interaction force during the motion 
of the observed particle. Thus, the interaction force en- 
tering Eq. ^ through becomes constant in this limit. 
The apparent entropy production then becomes equiv- 
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alent to that of an effective one-particle system subject 
to a Markovian dynamics with mean local velocity i>^, 
which trivially obeys the FT. This effective description 
is valid only to lowest order in time since taking into 
account higher order terms would include contributions 
arising from the correlated motion of the observed and 
the hidden particle. 

Although we are able to qualitatively understand the 
influence of coupling it remains a surprising feature why, 
in all of our experiments presented so far, only the slope 
of the FT is affected by the coupling strength while the 
linear relation Eq. ([T]) remains untouched. In order to 
elucidate this result we define the function 



/(AStot) = In [p(A5tot)/p(-AStot)] 



(7) 



which we assume to be analytic. First, we note that / 
is antisymmetric by construction, and thus for small en- 
tropy productions, Astot ^ 1, / trivially must be linear 
up to corrections of third order or higher [27| . Second, 
we discuss / for large entropy productions, Astot 3> 1. 
Solving Eq. ([7]) for p(— Astot) and integrating over all 
Astot yields 



-[-CX3 



p(AStot)e--^(^^~'-MAStot = 1, 



(8) 



by normalization. We assume that p(Astot) does not de- 
cay faster than a Gaussian as we have observed in all our 
measurements. For any quantity consisting of indepen- 
dent contributions the central limit theorem would imply 
a Gaussian. Any correlation will typically lead to an even 
slower decay. Under this assumption convergence of the 
integral in Eq. ^ requires that /(Astot) = 0(As'^^^). 
Since, in addition, / is antisymmetric we expect the 
generic asymptotic behavior to be linear, /(Astot) ^ 
Astot, with a slope generally different from the one for 
small Astot- 

Summarizing these arguments, we expect a linear func- 
tion both for small and for large entropy production for 
any time t. For intermediate entropy production this 
reasoning leaves the possibility of a nonlinear behavior. 
Even though we have found a constant slope for most 
experimental parameters, by fine-tuning the system and 
plasma parameter, we can observe an obviously nonlinear 
result, as shown in Fig. Sl^a). 

In contrast to the previous data, here, the two par- 
ticles are subjected to quite different potentials whereas 
the driving forces remain untouched. The potential of the 
hidden particle is adjusted such that it circulates freely 
along the torus {V2 = 71 keT) whereas a deep minimum 
{Vi = 262 ksT) remains in the tilted potential C/i of 
the observed particle, which, for F = 0, it is not able to 
leave. The latter's motion sets in only when the coupling 
helps it to surmount the potential barrier. This mode is 
identified in a typical trajectory shown in Fig. UJJb). The 
hidden particle (gray/thin line) moves with a period of 
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FIG. 4: (color online) (a) in [p(Astot)/p(— Astot)] as a func- 
tion of Astot for t = 3 s. (b) Typical trajectories of the ob- 
served (blue/thick line) and hidden particle (gray/thin line). 



4 s and almost constant velocity along U2{x2)- Around 
X2 = 0.75 it slightly slows down due to interaction. The 
reaction of the observed particle (blue/thick line) is more 
pronounced since locally it is confined within a potential 
minimum at Xi = —0.35 and there the interaction forces 
are dominant displacing it along positive xi. The appar- 
ent oscillations in the trajectory originate from the fact 
that not every time the observed particle is pushed (by 
the hidden one) this action results in a surmounting of 
the potential barrier. In most of the cases, the particle 
just relaxes to its original position. We observe that non- 
linearities in the intermediate regime of Astot are most 
pronounced when the trajectory length matches approx- 
imately the oscillation period. 

Concluding Perspectives. — We have investigated the 
influence of hidden slow degrees of freedom on the FT. 
In our experiments, we typically find that a FT-like sym- 
metry is preserved, however, with a different slope which 
depends, in particular, on the length of the observed tra- 
jectories. Consequently, in any experiment, where hidden 
slow degrees of freedom cannot be ruled out a priori, an 
observed linear behavior cannot be used to infer quanti- 
ties by implicitly assuming a — I. Theoretically, we have 
argued that a slope of 1 is to be expected for short tra- 
jectories while, for any length, for both small and large 
entropy production a constant slope should be generic. 
Classifying theoretically the conditions for finding an al- 
most constant slope over the full range, as we often did 
in our experiments, remains a task for future work. Like- 
wise, it will be important to explore, both in theory and 
experiments, how hidden slow degrees of freedom affect 
other quantities, like work and heat, their exact relations, 
and the fluctuation-dissipation theorem for NESS. 
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